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Abstract 

It is found that CP symmetry may be explicitly broken in the Higgs sector of 
a supersymmetric model with two extra neutral gauge bosons at the one-loop 
level. The phenomenology of the model, the Higgs sector in particular, is studied for 
a reasonable parameter space of the model, in the presence of explicit CP violation 
at the one-loop level. At least one of the neutral Higgs bosons of the model might 
be produced via the WW fusion process at the Large Hadron Collider. 
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I. Introduction 



Anticipations among high energy physicists for the discovery of new physics at the Large 
Hadron Colhder (LHC) are very high as it prepares to operate in full swing. There are a 
number of compelling rationales for anticipating new physics beyond the Standard Model 
(SM). One of them is the observed baryon asymmetry of the universe, which indicates the 
survival of more matter than antimater during the evolution of the universe. In the SM, 
the only source of CP violation is the complex phase in the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. It has been noticed that the size of CP violation in the SM by means of 
the CKM matrix alone is too small to explain the observed value of the baryon-to-entropy 
ratio, ub/s ~ 8 x 10~^^ [1], if the universe had begun from a baryon-symmetric state. 
Thus, in order to explain the observed baryon asymmetry of the universe, other sources 
of adequate CP violation are required. 

A number of alternative models beyond the SM have been investigated for the pos- 
sibility of CP violations. Supersymmetry (SUSY) has been with us for several decades, 
which nowadays is regarded as the most certain candidate for new physics. In reality, the 
necessity of CP violation beyond the SM is not the only raison d'etre for the SUSY. There 
are multiples of arguments that support its existence in nature. Some supersymmetric 
models have also been studied in this context, as their sophisticated Higgs sectors may 
possess sources of CP violation [2] . For some phenomenologically realistic supersymmetric 
models extended from the SM, soft SUSY breaking terms are essential ingredients [3]. If 
these soft SUSY breaking terms contain complex phases, the phenomenological analyses 
of these supersymmetric standard models might not only be complicated but also involve 
CP violation. 

The minimal supersymmetric standard model (MSSM) is the simplest version of su- 
persymmetric extension of the SM. Its Higgs sector has two Higgs doublets in order to 
give masses to up-like quarks and down-like quarks separately. At the one-loop level, a 
complex phase in the soft SUSY breaking terms of the MSSM induces an explicit CP 
mixing between scalar and pseudoscalar Higgs bosons [4]. 

Non-minimal versions of supersymmetric extension of the SM have additional Higgs 
singlets and thus can dynamically solve the dimensional /i-parameter problem in the 
MSSM by means of the vacuum expectation value (VEV) of the Higgs singlet [5,6,7]. 
They have also been studied within the context of explicit CP violation in their Higgs 
sectors [8,9,10,11]. The next-to-minimal supersymmetric standard model (NMSSM) is a 
typical member of them. Unlike the MSSM, the Higgs potential of the NMSSM has one 
nontrivial CP phase after redefining the Higgs fields at the tree level [9]. At the one- loop 
level, it also develops CP violating phases. The effects of explicit CP violation at the 
one-loop level in the NMSSM on the masses of neutral and charged Higgs bosons are 
predicted in the literature [10]. 

The Higgs potentials of both the minimal non-minimal supersymmetric model and the 
[/(l)-extended supersymmetric model may not have any CP phase at the tree level [11]. 
However, these models may also possess complex phases to induce explicit CP violation 
at the one-loop level, by taking the radiative corrections due to the quark and squark 
loops into account. 
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In this article, we would like to continue to study the possibility of CP violation in 
the Higgs sector of a supersymmetric Eq model. This model has two U{1) symmetries in 
addition to the SM gauge symmetry, thus with two additional neutral gauge bosons, and 
two Higgs singlets as well as two Higgs doublets [12,13]. The tree-level Higgs potential 
of this model may not have complex phase, because any complex phase can always be 
eliminated by rotating the relevant Higgs fields. At the one-loop level, it is shown that 
this model may allow CP violation in an explicit way due to radiative corrections. We 
study the Higgs phenomenology of this model by varying all the relevant parameters 
within reasonable ranges, to obtain the upper bound on the lightest neutral Higgs boson 
mass. We investigate prospects for discovering the neutral Higgs bosons of this model at 
the LHC, by calculating the minimum cross section for producing at least any one of the 
neutral Higgs bosons of this model via the WW fusion process at the LHC. 

II. Higgs Sector 

Let us describe the Higgs sector of our model. We assume that the electroweak gauge 
symmetry of our model is SU{2) x f/(l) x f/i(l) x U2{1), where the two extra U{1) 
symmetries are decomposed from Eq. Thus, it is a kind of rank-6 supersymmetric model. 
We assume that in general f/i(l) and U2{1) would mix with a certain mixing angle 6 to 
become two linearly orthogonal combinations, 11(1)' and U{1)". The Higgs sector of our 
model consists of two Higgs doublets, Hf = {Hl,Hi) and H2 = {H2,H2), and two 
neutral Higgs singlets, A^^i and A^2- The Yukawa interaction between Higgs superfields 
and quark superfields in the superpotential of our model may be expressed as [12,13] 

w ^ htQ^n2t% - hQ^nib% + min2Mi , (i) 

where we take only the third generation into account and ht and are respectively the 
dimensionless Yukawa coupling coefficients of top and bottom quarks, A is a dimensionless 
coefficient, Tii and 0,2 are the Higgs doublet superfields, N is the Higgs singlet superfield, 
and 6^ are respectively the right-handed top and bottom quark superfields, and Q is the 
left-handed SU{2) doublet quark superfield of the third generation. This superpotential 
has the same expression as discussed in Ref. [13] or Ref. [14], where relatively well-known 
rank-6 SUSY models are investigated. 

Note that, shown as the last term in the superpotential, only A^^i participates in cou- 
pling to the Higgs doublets, because the underlying Eq gauge symmetry does not allow the 
other Higgs singlet to do so [13]. Effectively, the coupling between A^i and the Higgs 
doublets corresponds to the fi term in the MSSM where the /i-parameter is generated by 
the VEV of the A^i. 

The Higgs potential of our model at the tree level is collected from D-terms, F-terms, 
and the soft terms in the superpotential. The most general form of the Higgs potential 
at the tree level is given as [13] 

Vo = mjHlHi + mlHlH2 + m^iV/iVi + m^iV|iV2 - {\AH1H2N1 + H.c.) 

+ \\\^[hIhihIh2 + hIhMni + HIH2NIN1] 
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lAI 



\HlH, 



|2 ^ dlA^HlH, 



HIH2) 



NlNi 



hNlN2) 



NlN2)f 

+ VlbCeiHlHi - nINi + N^N^)]^ , (2) 

where g2, gi, g'l, and g'^ are respectively the SU{2), U{1), U{1)', and U{1)" gauge couphng 
coefficients, A is a massive parameter, Cg = cos6 and Sg = sinO, and rrii {i = 1,2,3,4) 
are soft SUSY breaking masses. These four soft masses in the Higgs potential would 
eventually be eliminated by means of the minimum conditions for the Higgs potential 
with respect to four neutral Higgs fields. 

The parameters of the Higgs potential are assumed to be generally complex. Thus, A 
and A in the tree-level Higgs potential may be complex such that their complex phases 
may be factored out explicitly as AAe"^. We also assume that the VEVs, which four 
neutral components of the Higgs fields acquire after electroweak symmetry breaking, may 
in general be complex. However, by redefining the phases of Hi, H2, and N2, we may 



adjust the vacuum expectation values as vi = {Hf), V2 = {H2), Xie 



(iVi) and 



X2 = {N2), where vi, V2, xi, and X2 are real and the complex phase 0i is the overall phase 
in {H1H2N1) . Thus, looking at the Higgs potential at the tree level, one can easily notice 
that the only possible source of complex phases is XAH1H2N. By further redefining the 
phase of the Higgs singlet A^^i, it is always possible to make the two phases and 0i cancel 
each other so that the tree-level Higgs potential can be made completely real. Therefore, 
our model conserves the CP symmetry at the tree level. 

After the electroweak symmetry breaking, the tree-level mass of top quark is given as 
(^4^2)^, and the tree-level masses of stop quarks are given by the on-shell Lagrangian 



as 



m 



{ml + m%) 



+ mf: + -rriy cos 2/5 + G, =F \ Xt 



where mg and are the soft SUSY breaking masses for the stop quarks. 



VTir 



(3) 



gli)v^l2 with 
V2/V1, and 



VI + V2 is the squared mass of the neutral weak gauge boson, tan (3 



m 



T) 



2 2 



+ m;{Af + X-x\ coV(3 



2\AtXi cot P cos 




cos^ (3 - -VlOSgxi 



cos' p - -VlOCgxi 



4 




v"^ sin^ (3 




(4) 



with = (y'|f^/2 being the squared mass of the charged weak gauge boson, At being 
the trihnear soft SUSY breaking parameter of the stop quarks with mass dimension, and 
(j)t being a complex phase determined by 0i and the complex phase of At. Note that Gj is 
the effect of the two extra U{1) symmetries, but it does not contribute the mass splitting 
between the two stop quark masses. The mixing, and hence the mass splitting, between 
the stop quark masses is triggered by Xt. 

Now let us consider the one-loop radiative corrections to the tree-level Higgs potential. 
In supersymmetric models, the incomplete cancellation between ordinary particles and 
their superpartners yield the one-loop corrections to the tree-level Higgs boson masses. 
In SUSY models, the most dominant part of the one-loop corrections to the tree-level 
Higgs potential come primarily from the top and stop quark loops. For large tan/5 as 
large as 50, the contribution of the bottom and sbottom quark loops can also be large. In 
this paper, we consider the contributions from the top and stop quark loops at the one- 
loop level. The full Higgs potential at the one-loop level may be written as ^ = Vq + V^i, 
where Vi is contribution from the radiative corrections due to the top and stop quark 
loops. The effective potential method provides us [15] 



niMf 
647r2 



log 



A2 



(5) 



where A is the renormalization scale in the modified minimal subtraction scheme, the 
subscript I stands for the top and stop quarks: t, ii, M.i are the top and stop masses 
at the tree level given as functions of Higgs fields, and rii are the degrees of freedom for 
these particles. Including the sign convention, they are given as nt = —12 and rii^ = 6 
{i = 1,2), as in the above formula enter the stop quarks with a negative sign while the 
top quark with a positive sign. 

Since the parameters of the Higgs potential are assumed to be generally complex, we 
may have 0, which is the phase of XA. Unlike the tree-level case, we cannot redefine 
the phase of A^i at the one-loop level to cancel it. Thus, (j) may persist at the one-loop 
level. This can be most clearly be seen in the non-trivial tadpole minimum condition with 
respect to the pseudoscalar component of the Higgs field: 



= Asm, 



SmfAt sin i 



m 



t2' 



(6) 



IGvr^w^ sin^ 

where the first term comes from the tree-level Higgs potential and the second term comes 
from the radiative corrections, and the dimensionless function / arising from radiative 
corrections is defined as 

2 ™2" 



[mt 



-ml) 



m: 



m. 



log — - 



(7) 



But for the radiative corrections, the above tadpole minimum condition at the tree level 
would be satisfied when = 0. Due to the presence of the one-loop corrections, = is 
no longer in general the solution to the above tadpole minimum condition. 
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Our model has twelve real degrees of freedom in the Higgs sector. They may be 
classified as three neutral Goldstone bosons, a pair of charged Goldstone bosons, five 
neutral Higgs bosons and a pair of charged Higgs bosons. After the electroweak symmetry 
breaking, the three neutral Goldstone bosons and a pair of charged Goldstone bosons will 
be absorbed into the longitudinal component of Z, Z', Z" and W gauge bosons, where 
Z' and Z" are the extra neutral gauge bosons. 

The squared mass matrix M of the five neutral Higgs bosons is given as a symmetric 
5x5 matrix, obtained by the second derivatives of the Higgs potential with respect to 
the five neutral Higgs fields. At the tree level, the five neutral Higgs bosons may have 
definite CP parity, since the CP symmetry is conserved in the Higgs sector. Thus, we 
may denote them as Si (i = 1,2,3,4) for neutral scalar Higgs bosons and P for neutral 
pseudoscalar Higgs boson. In the {Si, S2, S3, S^, P) basis, the 5x5 matrix M at the tree 
level may be expressed as 



M = M^ + M' 



0' 



where M°' comes from the D-terms due to two extra U{1) symmetries of V'^, and M° 
comes from the remaining terms in Vq, namely, D-terms due to the SM gauge symmetry, 
the F-terms, and the soft terms of the tree-level Higgs potential. They may be expressed 



as 
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(9) 
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Explicitly, the matrix elements of MO and MO' are respectively given as follows: 



and 



MO, 
Mf 



M 



22 



33 
MO, 

MO3 
Ml, 

Ml, 



rr?2 cos^ /5 + Ml, sin^ cos^ a 
m^z sin^ /3 + Ml, cos^ (3 cos^ a 



m| /2) sin 2(3 - Ml, cos /? sin /3 cos^ a 



Ml, sin^ a 

2}?vxi cos [3 — Ml, sin [3 cos a sin « , 
2\^vxi sin [3 — Ml, cos 13 cos a sin a , 

2\Av—^ , 
sm 2a 



g'l'C-iy cos' (3 
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{g^ - g?)CeSev\os^ 13 



m: 



m: 



25 / '2^2 



^^^^hg',' - g?)CeSex\ 



9 

'^^Sg^Cl + g^Sl)x\ + \{g^Sl + g'-Cl)x\ + ^(^?? - 
18 o 9 



M- 



\QcI + g.'Siy sin 2/? - ^(^^l^ - g[')CeS,v^ sin 2/5 



^ / '2 ,^2 "2 rY2 \ n ^ / '2 rY2 

Y^l^i Cq + ^1 Sgjrai cos/? - -(^1 



+ 



{g\ - gi^)CgSevxi cos/3 , 



M' 



T^igiCg+gi SQ)vx2Cosp + -{g^SQ+ g^ Cq)vx2CosP 

18 D 

— {g'l - gi^)CgS0vx2 cos/3 , 



m: 



m: 



0' 
23 

0' 
24 



+ 



- y (fi'i^C'e + fi'?'S'g)rai sin/3 + '^—^{g'l - gl'^)CeSevXi sin/3 



- y + ^?'5g)ra2 sin/3 



2V15 



{g'l - g7)CeS0VX2 sin/3 
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(12) 



where tana = (?;/2xi) sin2/3 stands for the sphtting between an extra U{1) symmetry 
breaking scale and the electroweak scale. 

Note that both M° and M°' do not mix Si with P. In other words, there is no scalar- 
psuedoscalar mixing at the tree-level, hence the CP symmetry. It is straightforward to 
recognize that the single element M55 is the squared mass at the tree level of the neutral 
pseudoscalar Higgs boson. Note also that if the two extra U{1) symmetries are absent, 
we would have = 0. In this case, one of the neutral scalar Higgs bosons would be 
massless at the tree level, since M° may be decomposed into a block diagram consisting 
of three blocks, namely, 3x3 submatrix, M^, = and M55. 

Now, at the one-loop level, the squared mass matrix M of the five neutral Higgs bosons 
is corrected as 



M = M° + M°' + , 
where is the radiative corrections obtained from as 



(13) 
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(14) 
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Explicitly, the matrix elements of are given as follows, after imposing tadpole mini- 
mum conditions: 

, ^1 2 ■ 2 n 2 3 cos^ P ( Am 

M^^i = sm p cos a — 



Ml, 



Ml 




3 sin 2/? / Amy^ 5m| 



'55 A~ o n ■ n o (2 2 \2 1 

QTT^v^ sm p cos^ a [mf — m 



2 /loo 2 



= ~ ''^A "^os /3 sin /? cos^ a + 



327r2t;2 \^ 3 6 
3 sin/? frnf XxiAi^ cos/JA^-^ (ml At A 



vr^f^ \^ sin/3 2 y y sin^/3 

3 sin 2/5 4mj ^2 , 2\ f mfXxiA 



327r2t;2 \^sin2/3 ^ / V sin 2/3 

3 sin 2/3,^^ 2 2x [ml At A 



-{AGav' + m'z) 



3 sin 2/3,^^ 2 2\/ 4m? n 2\ 1 
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2 ■ n ■ SmfX'^Xi COS f3 2 2^ 



Svr^wtan/? \^ sin/3 ' 2 J (m?_^— m?^}^ 



327r2wtan/3 ' ^^(m|-m|) 



^vr^t; Y sin /3 ^ y (m?^ — m?^ ) 

3GeXiCOS/3 2^ 2^, 



Stt^i; y sin/5 2 y (m?_^— m?^) 

3GrfX2Cos/5 2^ 

3m^AMfxjA,-^sin(/>t ^(m?^, m|) 
Svr^f 2 sin^ /5 cos a (m?_^ — m?^ )2 
37nf\At cos sin (f)t g{m\, 

167r2t>tan/3sina {m~^ — m'^^y 
3m^\AtCospsmcl>t ^^^ ^ log(m|/m? ) 

327r"^t' tan p sm a (m|_^ — ; 



— m\ cos 13 cos a sin a 



SmjXAi^ fm'^AtAi^ ^sm/3Ai\ (7(m|, m|) 



Syr^ftan/^ \^ sin/3 ' 2 ^ (m?_^— m?^)^ 
_^ 3mgAcos/3A,-^ / 2mg m| _^ ^^^^2^ logl^/^^ 



lOvr^f lsin^/3 2 I {rrif — m? 



t2 tl' 

3^X1 f m^AtAj^ smpAi \ log(m|/m? ^ 
87r2t; \ sin/3 2 , 

3(j'c2;i sin d ( Arrii , ^ , 





sin/3AA log(m|/m? ) 



2 / 2 
{ml -ml] 



3GdX2smp Am^ ^,^2 2] ^ r^h'^ 



2 ^2 



+ " % — V;^ + 4Gbi;2 - m| log 



t2 



32;^V^ ""^^yl^^^V^ 
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3mf\ At COS sin (j)t gijn]^, mf^ ) 
IQir'^vsina (m? —rn^Y 
, 3mgA^cos/3sin0f Amj ^ ^ ^ \ log(m|/m|) 

M34 = o_..___ ^ „2\ + TT-. log' ' 



3^ SirHanp {ml - ml ) Stt^ ^\ 

SmfX'^AtXiAr sin 04 g{m1 , ml ^ 



Sn^v sin^ /3 tan P cos a (m| — m^J^ 
'''^iGcAAf cos^ /5 sin (f)t log(m|/m|~ 



+ 



STT^tanasina (m? — m? ) 



,.1 _ Sm^GdAtX2\ sin 0^ log(m,-^/mgJ 



Svr^ tan sin a (m?__^ — m?^ 



where 



2 SAm^AtCos^t 2 2N 

= ^ 9-7:/ , fnt ) , 16 

^ 87r2t;sin2asin2/3'' ^ ' ^ ^ 

Af^ = Aa; cot P — At cos 0t , 

= Ax cot P cos </)t — y4t , 

A, = f _ - ,4) + f l!^ _ cos 2,4 . (17) 



G. = |^(4C2» - 1) - |J(4Cm + 1) , 

'2 "2 

= %:{VlES2e + C2e-4)-^{VlES2e + C2e + ^) , 
6b 6b 

'2 "2 

io io 

'2 "2 

Grf = 1^(10 -3v^S2, + 5C2e) + 1^(10 + 3v^52e - 5C2,) , (18) 
and the dimensionless function g is defined as 

^?K,rnJ) = ^±^log^ + 2. (19) 
" m^ — my m^ 

Note first that the matrix elements M^g {i = 1,2,3,4) are proportional to sin0j. If 
0f = 0, these elements would be zero, and the scalar-psuedoscalar mixing at the one- loop 
level would not occur in the Higgs sector. Therefore, there would be no CP violation in 
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the Higgs sector at the one-loop level. The squared mass of the pseudoscalar Higgs boson 
at the one-loop level, mp, would be given simply by the (5, 5)-th element of the A-f, taking 
0i = 0. It is given by adding the radiative corrections as 

ml = 2XAv 9 ^ / , "^f ■ 20 

^ sin2a 87r2t; sin 2a sin^ /^'^ ^ ^ ' 

In this case, the D-terms of extra U{1) symmetries would not contribute to the mass of 
the pseudoscalar Higgs boson either at the tree level or at the one-loop level. 

If, on the other hand, (pt 7^ 0, there would be CP violation at the one-loop level. The 
CP phase in the radiative corrections generates the scalar-pseudoscalar mixing, thus the 
five neutral Higgs bosons are no longer states of definite CP parity. In this case, the 
mass matrix should be diagonalized to obtain mass eigenstates hi {i = 1, 2, 3, 4, 5) whose 
squared masses m|, (i = 1,2,3,4,5) are the eigenvalues of the mass matrix. These five 
neutral Higgs bosons are usually numbered such that hi is the lightest neutral Higgs boson 
and /is is the heaviest. Hereafter, we work in the explicit CP violation scenario, that is, 
with 04 7^ 0. 

In our model, the squared masses of the two extra gauge bosons m^, and Tfi ^11 are 
obtained as the eigenvalues of the mass matrix for them. The explicit expressions for m^; 
and m^// are given as 

1 



^l" = + m|J + ^m%^ - m|)2 + 4A2 , (21) 

where 

= ^g[\\4-C2e + Vl5cos2pS2e) + ^g?xlS's 

+ ^g'M{S + 7C2e-Vl5S2e) , 
3d 

3o 



Ai = ^g[mzv{\/5 cos2f] Se + V^Ce) 
= -gimzv{V5cos2(3C0 - VsSe) 



:22) 



The two mixing angles in our model, ai between Z and Z' and 02 between Z and Z", are 
expressed as 



1 _i / 2A, 

5 



ai = - tan — , (23) 

2 V mi, — mi / 



for i = l,2. 
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III. NUMERICAL ANALYSIS 



For the sake of simplicity, we take g'^ = g[ = ^JbJSgi in our numerical analysis, motivated 
by the gauge coupling unification. We consider the region of the parameter space bounded 
as < 61 < 7r/2, < 0t < tt, 1 < tan/? < 30, and < A < 0.83. We assume that the 
lighter stop quark is heavier than the top quark. We also assume that all of the relevant 
mass parameters, mp, mq, rriT, and At, vary within the range of 100 to 1000 GeV. Note 
that we employ the mass of the pseudoscalar Higgs boson at the one-loop level in the CP 
conserving scenario, mp, instead of A as an input parameter. Further, we use a combined 
constraint of Axi > 150 GeV, as the experimental data on the chargino system set the 
lower bound on the effective fi parameter, fi = \xi. For the values of xi and X2, we would 
set their ranges not by hand but by experimental constraints. 

There are strong experimental constraints on the mass of the extra neutral gauge 
boson and the mixing between Z in the SM and the extra neutral gauge boson. Thus, 
any model with extra neutral gauge bosons, such as our model, should comply with these 
constraints, whose exact values may dependent on the specific structures of the models. 
We would like to take in this article that the mixing angles, ai and a2, should be smaller 
than 3 x 10""^ and the masses of the two extra gauge bosons, m^/ and m^//, should be 
larger than 800 GeV. 

The experimental constraints on the Higgs sector should also be taken into account 
in the numerical analysis. The latest experimental analyses tell that the SM Higgs boson 
lighter than 114.5 GeV is excluded at the 95 % confidence level. This lower bound on 
the SM Higgs boson mass may be applied to our model by considering the relevant Higgs 
couplings. Recently, the LEP collaborations reported the model-independent upper bound 
on [qzzh / Q^hY 9^ % confidence level [16]. 

First, we determine the values of Xi and X2, varying the values of other relevant pa- 
rameters within their allowed ranges. While the values of the above parameters are chosen 
by the random number generation method within their respective ranges, the values of Xi 
and X2 are determined in terms of the other parameters by imposing experimental con- 
straints. The result is shown in Fig. 1(a), where a distribution of 7125 points is displayed 
in the (a:i,a:2)-plane. These points are selected among 10^ random points as they satisfy 
all of the above experimental constraints. Each point represents a set of parameter values, 
of which the values of Xi and X2 are explicit while others are implicit. It is notable that 
these selected points are distributed in the area of the (xi,X2)-plane where X1 + X2 > 2100 
GeV. Some points are scattered at xi ~ 400 GeV. 

Then, we calculate for each point in Fig. 1(a) the mass of the lightest neutral Higgs 
boson in our model. In this way, it is clear that the results are consistent with the relevant 
parameter ranges as well as the experimental constraints. The result is shown in Fig. 1(b). 
It is quite remarkable that the majority of the points are scattered within the range of 
117 < nih^ < 140 GeV, while a few of them are distributed where m/ij is as low as 30 
GeV. The result of Fig. 1(b) suggests that the mass of the lightest neutral Higgs boson in 
our model is most probably about 130 GeV at the one-loop level. 

One may notice some pattern in Fig. 1(b). We find that this pattern comes from the 
experimental constraints on the extra gauge bosons rather than that the experimental 
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bound on the SM Higgs boson mass. The lower bound on Xi + X2 is found to arise from 
the experimental constraints on the masses of the extra neutral gauge bosons. Meanwhile, 
most of points with rrih^ < 115 GeV are excluded by the experimental constraints on the 
SM Higgs boson mass. We also calculate the masses of other neutral Higgs bosons. The 
results are shown in Figs. 1(c) and (d), where we display the correlation between m/13 
and m/i2 in Fig. 1(c) and the correlation between rrih^ and nih^ in Fig. 1(d). The 
points in these figures are obtained with the same parameter values as in Figs. 1(a) or 
(b). Note the clear hierarchy between the masses of the neutral Higgs bosons such that 
m^g > rrLfi^ in Fig. 1(c) and m/^g > rrih^ in Fig. 1(d). The ranges for the masses of heavier 
neutral Higgs bosons in our model, estimated using the aforementioned parameter values, 
are: 100 < ruh, < 997 GeV, 116 < rrih, < 998 GeV, 262 < 711^, < 1189 GeV, and 
987 < m/ig < 1536 GeV, where the upper bounds come from theoretical arguments and 
the lower bounds come from phenomenological constraints. 

Now, we examine the possibility of discovering one of the neutral Higgs bosons in our 
model in the pp collisions at the LHC, where the most dominant process for the Higgs 
production is the gluon fusion process, with thick QCD backgrounds. The WW fusion 
process is considered as the next dominant process for the Higgs production, which is 
relatively cleaner than the gluon fusion process. We would like to focus on the WW 
fusion process. 

We find that the PYTHIA program is useful for calculating the Higgs production 
mechanism than for other processes, although it has not yet been applied to the CP 
violation scenario in the MSSM Higgs sector. However, the production cross section of 
the neutral Higgs bosons in our model with explicit CP violation via the WW fusion 
process in pp collisions is obtained by using the PYTHIA 6.4 program after appropriately 
modifying the relevant Higgs coupling coefficients [17]. More precisely, we normalize 
GwwH,, the WW hi coupling coefficient of the Higgs coupling to a pair of W bosons, by 
the corresponding SM Higgs coupling coefficient. We have 

Gwwh, = {cospOu + sm(302i) , (24) 

where Oij = 1, 2, 3, 4, 5) are the elements of the orthogonal matrix that diagonalizes 
the mass matrix for the five neutral Higgs bosons. 

Technically, we set the number of events to generate for each point as NEV = 2000. 
The Higgs coupling coefficient is set by MSTP(4)=1, and the normalized Higgs coupling 
to a W boson pair is set by PARU(165) = Gwwhi- The factorization scale and the 
renormalization scale are taken to be the neutral Higgs boson mass, that is, PARP(193) 
= PMAS(25,1) and PARP(194)= PMAS(25,1). The PDF library of the CTEQ5L is used 
in our program, MSTP(51)=7, which is the default parton distribution function set for 
the proton in PYTHIA 6.4. We use MSTP(33)=0 to include the K factor in hard cross 
sections for parton interactions in PYTHIA 6.4 by default. The WW fusion process for 
the lightest neutral Higgs boson is set by MSUB(124) = 1. 

In this way, we obtain all of (Twwhi {i = 1,2,3,4,5), the production cross sections of 
hi in our model with explicit CP violation via the WW fusion process in pp collisions. 
They are given as functions of the participating neutral Higgs boson masses. Among the 
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five production cross sections, we select the largest one, as we are interested in discovering 
any one of tlie five neutral Higgs bosons. Thus, we introduce 

<^WWh = ^-^'^{<7WWhi,<ywWh2,<^WWh3,0'wWh4,0-wWhs) ■ (25) 

We show our result in Fig. 2, where we plot (T\YWh against m/^^. The parameter values 
for each point are the same as in Fig. 1(a) or Fig. 1(b). We find that the smallest value 
for awwh is about 1 pb. This implies that at least one of the five neutral Higgs bosons 
in our model may be produced with its cross section larger than 1 pb. The accumulated 
integrated luminosity of 30~^ fb at the LHC would yield 6000 raw Higgs events, if we 
allow 20 % for the efficiency and acceptance. Therefore, we expect with relatively strong 
confidence that at least one of the five neutral Higgs bosons in our model might be 
produced via the WW fusion process at the LHC, if they exist. 

Here, the roles that the exotic quarks take part in are worth mentioning with respect to 
the Higgs phenomenology of our model. The exotic quarks may inhabit the fundamental 
27 representation of Eq, which is the underlying gauge symmetry of our model. In the 
fundamental 27 representation, 15 components are occupied by the SM matter fields, 
4 components by the two Higgs doublets, 2 components by the Higgs singlet, and the 
remaining 6 components are occupied by the exotic quarks [12-14,18-20]. 

The form of the superpotential tells us that the exotic quarks may couple to various 
Higgs fields. They couple directly to the neutral component of the Higgs singlet A'^i and 
indirectly, through the mixing among the neutral Higgs bosons via the diagonalization 
matrix, to other neutral Higgs fields. If the masses of the exotic quarks are comparable 
to the electroweak symmetry breaking or SUSY breaking scales, the low energy SUSY 
phenomenology might be affected by their presence. The effects of the exotic quarks 
might appear in the gluon fusion processes for Higgs productions, as well as in the Higgs 
decay processes. In particular, for example, the Higgs decays into a pair of gluons or 
photons might receive the effects of the exotic quarks, when the mass of the Higgs boson 
is below the range where the decay channel into a pair of gauge bosons are not yet open. 

However, it is somewhat difficult to predict the amount of the exotic quark effects be- 
cause it depends on the relevant parameters in a complicated way. The coupling strength 
of the exotic quarks to the neutral Higgs bosons are weak in general but might be strong, 
depending on what is the explicit structure of the orthogonal matrix that diagonalizes 
the mass matrix for the neutral Higgs bosons. Therefore, it would be valuable to study 
elsewhere a comprehensive research on the effects of the exotic quarks in our model. 

IV. CONCLUSIONS 

We study a supersymmetric Eq model with two extra U{1) symmetries besides the SM 
gauge symmetry, and two neutral Higgs singlets besides two MSSM Higgs doublets. We 
find that the Higgs sector of our model may generally accommodate a non-trivial complex 
phase which can cause the scalar-pseudoscalar mixing among the five neutral Higgs bosons, 
by virtue of radiative corrections due to the top and stop quark loops. Thus, explicit CP 
violation at the one-loop level is viable in our model. 
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Numerical analysis shows that there are parameter regions in our model which comply 
with a number of experimental constraints such as the lower bound on the extra neutral 
gauge boson masses and the upper bound on the mixing between the extra neutral gauge 
bosons and the SM neutral gauge boson. Within the allowed parameter regions, we study 
the behavior of the vacuum expectation values of the two Higgs singlets, xi and X2- We 
find that they cannot be simultaneously small. The experimental constraints on the extra 
neutral gauge bosons restrict that xi + X2 should be larger than 2100 GeV whereas either 
one of them may be as small as 400 GeV. 

The possibility of discovering one of the five neutral Higgs bosons in our model is 
examined by calculating the production cross sections using the PYTHIA 6.4 program, 
where the relevant Higgs couplings are modified suitably. We focus the WW fusion process 
at the LHC for their productions. We find that at least one of five neutral Higgs bosons 
can be produced enough via the WW fusion process at the LHC. Thus, we speculate that 
the present SUSY Eq model can be tested by the Higgs searches at the LHC. 
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FIGURE CAPTION 



FIG. 1(a). : A distribution of 7125 points in the (xi, X2) plane. Each point represents a set 
of parameter values that satisfies the experimental constraints on the extra neutral gauge 
boson masses, on their mixings with the SM neutral gauge bosson, and on the SM Higgs 
boson mass. The values of xi and X2 are explicitly shown, and the other parameters have 
certain values within their ranges respectively by the random number generation method: 
1 < tan/3 < 30, < A < 0.83, < ^ < 7r/2, < 0* < vr, 100 < mA,mQ,mT,At < 1000 
GeV. 

FIG. 1(b). : The plot of m?,^ against (xi + X2). For each of the 7125 points in Fig. 1(a), 
the mass of the lightest neutral Higgs boson is calculated in terms of the parameter values 
represented by the point. 

FIG. 1(c). : The distribution of 7125 points in the (m/jg, m/j2)-plane. They are distributed 
between 100 < nih^ < 997 GeV and 116 < mh,^ < 998 GeV, and they satisfy > m^^. 
These points are obtained with the same parameter values as in Figs. 1(a) or (b). 

FIG. 1(d). : The distribution of 7125 points in the (m/i^, m/i^)-plane. They are distributed 
between 262 < m/j^ < 1189 GeV and 987 < m^g < 1536 GeV, and they satisfy m^jg > m/^^. 
These points are obtained with the same parameter values as in Figs. 1(a) or (b). 

FIG. 2. : The polt of awwh against m/^^. The production cross sections of the five 
neutral Higgs bosons via WW fusion process in pp collisions are calculated in terms of 
the parameter values represented by the point, and the largest of them is chosen, for each 
of the 7125 points in Fig. 1(a). 
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FIG. 1(a): A distribution of 7125 points in the (xi, X2) plane. Each point represents a set 
of parameter values that satisfies the experimental constraints on the extra neutral gauge 
boson masses, on their mixings with the SM neutral gauge bosson, and on the SM Higgs 
boson mass. The values of Xi and X2 are explicitly shown, and the other parameters have 
certain values within their ranges respectively by the random number generation method: 
1 < tan/3 < 30, < A < 0.83, < ^ < 7r/2, < 0t < vr, 100 < mA,mQ,mT,At < 1000 
GeV. 
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FIG. 1(b): The plot of m/i^ against {xi + 2:2). For each of the 7125 points in Fig. 
the mass of the hghtest neutral Higgs boson is calculated in terms of the parameter 
represented by the point. 
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FIG. 1(c): The distribution of 7125 points in the (m/jg, m/ij) -plane. They are distributed 
between 100 < nih^ < 997 GeV and 116 < nih-^ < 998 GeV, and they satisfy > mh^. 
These points are obtained with the same parameter values as in Figs. 1(a) or (b). 
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FIG. 1(d): The distribution of 7125 points in the (m/i^, m/j^) -plane. They are distributed 
between 262 < rrih^ < 1189 GeV and 987 < rrih^ < 1536 GeV, and they satisfy ruh^ > rrih^. 
These points are obtained with the same parameter values as in Figs. 1(a) or (b). 
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FIG. 2: The polt of awwh against m/^^. The production cross sections of the five neu- 
tral Higgs bosons via WW fusion process in pp coUisions are calculated in terms of the 
parameter values represented by the point, and the largest of them is chosen, for each of 
the 7125 points in Fig. 1(a). 
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